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£SJ \ We consider the Maxwell-London electrodynamics of three dimensional superconductors in p- 

» I . wave pairing states with nodal points or lines in the energy gap. The current-velocity relation is 

then nonlinear in the applied field, cubic for point nodes and quadratic for lines. We obtain explicit 
angular and depth dependent expressions for measurable quantities such as the transverse magnetic 
moment, and associated torque. These dependences are different for point and line nodes and can 
be used to distinguish between different order parameters. We discuss the experimental feasibility 
of this method, and bring forth its advantages, as well as limitations that might be present. 
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I. INTRODUCTION 
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^ The number and variety of superconducting materials for which evidence of exotic Cooper pairing (i.e., pairing in 

a state other than the usual s-wave) exists is constantly increasing. For high temperature superconducting oxides 
(HTSC's) the consensus!!! is indeed that the pairing state is, in nearly all cases, at least predominantly d-wave, 
specifically of the d x 2_ y 2 form, with lines of nodes- Rather persuasive (although not conclusive) evidence in the form 
\ of both experiments and theoretical arguments J3~u has recently been brought fa-ward for p-wave superconductivity 
in Sr2Ru04. The pairing state currently favored by many is of the same formETQ as that of the A phasda in 3 He, 
which has point nodes. Several heavy fermion (HF) materials, the discovery of which predated that of HTSC's but 
for which determinations of the pairing state have, proved harder to achieve, are now also, believed with varying 
degrees of certainty, to belong in the exotic campE£nl3. There are also results indicatingoo that superconducting 
t-H ! families of organic salts such as k - (BEDT - TTF) 2 Cu(NCS) 2 and (TMTSF) 2 X (X = PF 6 , C10 4 , etc.) also e: 
^ ■ unconventional superconductivity. In some cases it has been argued that the pairing appears to be in the p-wavi 

Determination of pairing states is not easy, particularly if one wishes to know more details than merely their 
overall symmetry. Even in the best studied HTSC's, questions such as what is the angle between lines of nodes in 
orthorhombic compounds, or whether true nodes, rather than very deep minima, exist, are still matters for occasionally 
heated debate. The situation is much pppp f° r the other materials mentioned, where the evidence is much more 
preliminary, and at times contradictory.EZrE3 The determination of the pairing state is often hampered by difficulties 
in interpreting results. Regions (points or lines) where the energy gap vanishes are often the signature of exotic pairing 
(but not invariably, the B phase of 3 He is a well-known counterexample). These "gap nodes" lead to various power 
law behaviors for quantities that otherwise would behave exponentially with temperature, but sometimes there are 
alternative explanations for the power laws. It is difficult moreover, to distinguish between experimental outcomes 
arising from zeroes in the energy gap and those arising only from strong anisotropy. An additional complication is 
that for non s-wave superconducting materials, the order parameter (OP) state at the surface may easily differEj from 
q ' that in the bulk. 

(j , It is therefore important to study probes of the OP symmetry able to discern as unambiguously as possible details 
of the pairing state, such as the existence, nature and position of the bulk OP nodes. One such probe is afforded by the 
nonlinear Maxwell-London electrodynamics of exotic pairing states in the Meissner regime. Electrodynamic effects 
probe the sample awer a scale determined by the penetration depth A, which is large for the materials of interest. 
It was pointed outE-3 in the context of d-wave superconductivity that order parameter nodes lead to observable 
nonlinear effects at low temperatures, the chief quantities of experimental interest being the magnetic field dependent 
penetration depth X(H), the nonlinear transverse component of the magnetic mnmcnt, mi, induced by the application 
of a magnetic field, and the torque associated with this transverse moment riTrJI Further developments of the method, 
always in the context of predominantly d-wave superconductivity, showedEa that it can be used to perform node 
spectroscopy, that is, to infer in detail the angular structure of the regions where the order parameter vanishes (nodes) 
or is very small ("quasinodes"). 

These developments took place within the study of the high temperature oxide superconductors. For these materials, 
the temperature scales as set by T c are higher and achieving the required low temperature conditions is very easy™ 
However, recent improvements in experimental techniques involving torsional oscillatorsEj and torque magnetometryE2l 
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make it possible to measure extremely small moments and torques at dilution refrigerator temperatures. Experiments 
to accurately measure A in that temperature range are alscEa being planned. Thus, the relevant region for performing 
nonlinear electrodynamics experiments in low T c materials is becoming accessible. 

With this in mind, we take up in this work the question of the use of methods based on nonlinear electrodynamics 
to study exotic superconducting materials, other than HTSC's. Specifically, we will consider here simple OP's both 
with point nodes and with three-dimensional nodal lines, as would occur for example in p-wave superconductivity. 
Our efforts will focus on the calculation of the dependence of m± (or its associated torque) on the magnetic field and 
the appropriate angle of rotation. We also compute the field dependence of the low temperature penetration depth. 
We will present estimates based on published values of the relevant material parameters showing that the required 
measurements appear to be technically feasible. These estimates are presented, not to prejudge the pairing state 
associated with any material, but rather to show the expected signal if the material indeed does have the assumed 
OP. 

In the next Section we introduce the geometries and the order parameter forms that we study. We then calculate 
the nonlinear relation between current and superfluid flow field, using an extension of the three dimensional methods 
of Ref. E^. pErom these relations, we obtain the physical quantities of interest, through the appropriate generalization 



of existing^ perturbation methods. In Section HI we summarize our results, and consider the question of the 



experimental feasibility of using this method on several materials. We conclude with a discussion of the advantages 
and limitations of the method and of the specific treatment presented in this work. 



II. METHODS AND RESULTS 



A. Maxwell-London Electrodynamics 

We first briefly outline the nonlinear Maxwell-London equations, on which our method is built. We will not dwell 
into any details that were discussed elsewhere. ot3c3. When a magnetic field H a is applied to a supercoaductor a 
superfluid flow field v(r) is set up. The relation between v(r) and the local magnetic field H(r) is giver£3 by the 
second London equation: 

Vxv=-H. (2.1) 



where e is the proton charge. Ampere's law for steady-state currents, V x H = — j, can be combined with Eq. (2.1) 
to obtain: 

VxVxv=^j(v). (2.2) 

In this equation there are still two unknnwp fields. For a solution to be obtained, the functional relationship 
between j and v is needed. This can be founded by using the two-fluid model. The quasiparticle excitation spectrum, 
E(e) = (e 2 + lA(s)! 2 ) 1 / 2 , is modified by a Doppler shift to E(e) + v/ ■ v. Here e is the quasiparticle energy referred 
to the Fermi surface, A(s) denotes the OP dependence on the point s on the Fermi surfaceo, and Vf is the Fermi 
velocity. This leads to a relation between j and v of the form 

j(v)=j Kn (v)+j„i(v). (2.3) 
After some algebrariE, the linear and nonlinear parts can be written, respectively, as: 

jim(v) = -eN f / d 2 s n(s)v/(v/ ■ v), (2.4a) 

JFS 



j n/ (v) = -2eN f / d 2 s n(s)v f / de f(E(e) + V/ • v) 



FS 



(2.4b) 



where Nf is the total density of states at the Fermi level, n(s) the local density of stat es at the Fermi surface (FS), 
normalized to unity, and / the Fermi function. The first term in (2.2), given by ( [2.4a| ), is the usual linear r elatio n 
hin = — epv, where p is the superfluid density tensor. At T = the nonlinear (in v) corrections, described by (2.4h), 
can be written as: 
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j„Kv) = -2eN f / d 2 s n(s)v f Q(-v f ■ v - |A( S )|)[( V/ • v) 2 - |A( S )| 2 ] 1/2 , (2.5) 
Jfs 



which is valid for any A(s). The key point is that the step function in Eq.(2.5) restricts the integration over the FS 
by 

|A( s )|+ V/ -v<0. (2.6) 

Thus, when the OP has nodes (or very deep minima) in the Fermi surface, only regions near these nodes participate 
in populatingj-the quasiparticle spectrum. The integration is dominated by contributions from these regions, and can 
be written asE2l a sum over local contributions from each of them. The values of the Fermi velocity in the integrand 
can be replaced by their local values at the corresponding nodal region. We need, therefore, more information on the 
geometry of the superconductor and the angular dependence of the order parameter to carry out the above integration. 



B. Geometry and order parameters 

We will consider superconducting materials of orthorhombic or higher symmetry and denote the crystallographic 
axes as a, b and c, corresponding as customary to the x, y and z directions. We will assume that the samples are 
infinite in a plane parallel to the direction of the applied field, and of thickness d in the direction normal to this plane. 



This allows us to solve (2.2) and (2.5) analytically. Effects of the sample finite extension in the plane would have,-to 
be taken into account numerically, but this is unnecessary since it has been shown in the context of d-wave pairingu, 
that such effects merely lead to a small increase in the amplitude of the nonlinear signal, and to no change in its 
angular or field dependence. On the other hand, the effects of the thickness d are very important and we will include 
them fully. 

We will consider two simple tyrocp,of p-wave OP's in this paper. The first is representative of the case where the OP 
has point nodes, as might occumoQ in Sr2Ru04. Up to a phase factor (the nonlinear electrodynamic effects depend 
only on the absolute value of the OP) we write for the angular dependence of the OP near the nodes: 

A(0) = A o sin(0), (2.7) 



where 9 is a polar angle. Only the local properties at the nodes are important: the form (2.7) is assumed only near 
the nodes, e.g. near 6 = Q it means A(0) « A o 0. Thus, the parameter A must be thought of as the slope of the 
OP near the node, rather than its maximum value. The second-type we wilLconsider is a prototype of OP's with line 
nodes, as they might occur in some heavy fermion compoundsEj or even inEj Sr2Ru04. Again, up to an unimportant 
phase factor, we assume the form for the angular dependence near the nodal line: 

A(0) = A o cos(0), (2.8) 

where the above warning as to the interpretation of Ao as the slope near the nodal region must be repeated. These 
two forms are archetypes for the possibly more intricate forms of the angular dependence of the OP in real materials. 
In non s-wave superconductors., the OP need not belong to a one-dimensional representation. Because of this, there 
may be OP collective modesQ'Ea and internal structure effects that are not included in our considerations. The angular 
dependence of the OP, however, will be in a solid very strongly pinnecfl by crystal effects (this is obviously not the 
case in liquid 3 He). The internal structure of the OP should then not affect our nonlinear results, since only the 
application of small dc fields is involved. 

The experimental setup we envision would involve applying a field parallel to the a — c plane, with the direction 
normal to the slab being along the b axisEJ. The sample would then be rotated about the b axis while the magnetic 
field remains fixed. The currents then flow in various directions depending on their orientation relative to the OP 
under consideration. We will denote by ip the angle between the applied field H a and the z axis and we will investigate 
the angular dependence of the transverse magnetic moment or the torque as a function of tp. We will also calculate 
the field dependence of the penetration depth for the directions of symmetry. 

In this geometry we can solve the problem analytically. When H a is applied in the a — c plane, the fields have only 



x and z components, which depend only on the coordinate y. Eq. (2.2) then reduces to 



d 2 v Aire . . . . . 

+ — J v = 0. (2.9) 
ay z c z 

For our given geometry, j„h and Vi have odd parity with respect to the y coordinate, so it is sufficient to solve the 
boundary value problem for y > 0. The two required boundary conditions are: 
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- (Vxv)| y=d/2 =H a , 
v| a=0 = 0. 

We can now proceed to explicitly calculating the nonlinear currents. 



(2.10a) 
(2.10b) 



C. Nonlinear currents 



First, we carry out the integration in Eq. (2.5). This can be performed exactly on a three dimensional Fermi 
surface, without the need to take recourse to the approximations discussed in Appendix A of Ref. ^7]. The relevant 
regions of integration as discussed below (2.6) are contained within a small range near the nodes, with boundaries 
that can be expressed in terms of lim itin g angles 9 C , as determined from (v/ • v) 2 = |A(6> C )| 2 . 

Consider first the OP given in Eq.(2.7). In this case v/(s) can be replaced in the relevant regions of the integrand 
by its local value, along the z axis, at the nodes. By symmetry, we can restrict ourselves to the nod e at 9 = since 
the contribution from 9 = it is identical. Thus, we have v/ sa (0,0, Vf z ), and the restriction (2.6) means that in 



performing the integral in (2.5), we can replace J FS d 2 s n(s) by J n d<fi9d9 /4-7T, where f2 c denotes the region \9\ < 9 C 



with 9 2 = (vf z v z ) 2 /AQ, and with no restrictions on <fi. It is easy to see that this yields only a z-component to the 
nonlinear current. The integrals are elementary and one finds: 



Jnlz 



1±. 

A 2 



(2.11) 



where Aq the local gap slope. 



For the case of an order parameter as given in (2 



where the nodal line is at 6 = 7r/2, we can take Vf z = over 
the region of integration, which is then limited to \9 — tt/2\ < 9 C , where (9 C — it/2) 2 = (cosavfv± /A ) . Here v± is 
the projection of v on the x — y plane, and a the angle between v± and the in-plane vj. We make the replacement 
J FS d 2 s n(s) — > j n d<ftd9/4ir, where Q c is the region of integration as defined by 9 C . For orthorhombic symmetry, 
Vf depends on (f>. We then transform the integral over <p to one over a usin g the relation <j) = j3 + a, where (3 is the 
(fixed) angle v± makes with the x axis, and a is restricted (from Eq. (2J3)) to tt/2 < a < 3tt/2. The integration is 
lengthier but straightforward and results in two components to the current. The x component is: 



Jnlx 



1 AT "fx 



>V2 



Ax(v), 



(2.12) 



where A^v) = (l/5)[(3 + 2d 2 ) + (1 — S 2 )(v 2 /(v 2 + v 2 ))} (A x = 1 when the material has tetragonal symmetry), Vf x 
is the Fermi speed along the x axis and the a — b plane anisotropy is characterized by 8 = X x /\, where Xi denotes 
the pene tration depth along th e i d irectio n. Th e y component, j n iy, is obtained by making the obvious replacements 
in (2.12). The coefficients in (2.11) and ( 2.1 2| ) are given in terms of the local values of the Fermi velocity and are 
therefore independent of the detailed shape of the Fermi surface. 



These expressions for j as a function of v can be inserted into Eq. (2.9), which then becomes a nonlinear differential 
equation in terms of the flow field only. Implementing a perturbation scheme to lowest order in the flow field will 
allow exact expressions for both OP's to be obtained. This is addressed in the next subsection. 



D. Perturbation solution 



For an OP of the form (2/7), we can insert Eq. ( 2.11 ) into Eq. (2J3). We can then write the equation for the 
component carrying the nonlinear term as: 



dY? 



— v z = 0, 

Vr 



(2.13) 

where we have introduced the dimensionless coordinate Yi = y/Xi- The local critical velocity is defined as v C i = An /ffj, 
(for i = x,y,z) and we have used the three dimensional relation 1/A 2 = (47re 2 /3c 2 )N fv 2 ^ From Eq. ( 2.10a ), the 
boundary condition at the surface Yi = Yi S = d/(2Xi) can be written in terms of our new variables: 



dv x 
dY x 



eX x 



H a cosip, 



Y=Y X 



dv z 
dY z 



eA, 



-H a smip. 



(2.14) 



=Y za 
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We now expand v z (Y z ) to first order in the parameter a z = (f/ z /Ao) 2 , which is small in the typical experimental 
situations. We write v z (Y z ) = vq z (Y z ) + a z v\ z {Y z ). To zeroth order, we have the usual linear equation: 



d 2 v 0z 
dY? 



v 0z = 0, 



with vq z satisfying the boundary conditions (2.14), (2.10a). The solution is: 

v z(Y z ) = c z smh(F z ), 



where 



The nonlinear part v\ z satisfies: 



The boundary conditions are dviz/dY^ 



e\ z H a smijj 
ccosh(Y zs ) ' 



d 2 vi z 3 

- Vlz+Voz = . 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



elementary methods and is given by 

vu(Y z ) = {1/8)4 fa smh{Y z ) + 3Y Z cosh(F z ) - (1/4) sinh(3r 2 )] , 



= and vi z (0) = 0. The complete solution to Eq. (2.18) is found by 

(2.19) 



where c x = (3/2)(sinh(2Y; s ) - 2Y ZS ) tanh(y zs ) - 9/4. 

The magnetic field in the sample can be calculated from the field v via Eq. (2.1). Including also the purely linear 
component arising from v x we obtain: 



H X {Y Z 



H a sin V> 
cosh(y 2S ) 



cosh(y z ) 



1 / H a sin tp 
8 \H Qz cosh(Y zs ) 



Jh (Y z ) 



(2.20a) 



H Z {Y X ) = g "Jvt cosh(y x ), (2.20b) 
cosh(Yx S J 

where the nonlinear depth dependence is contained in 

MY,) = 3Y Z sinh(r z ) - (3/4) cosh(3F z ) + ( Cl + 3) cosh(y z ), (2.20c) 
and we have introduced the characteristic field 

m. - J^- (2.2D 

Here 4>o is the superconducting flux quantum, and & = u/j/(7rAo) is the local coherence length. As opposed to the 
d-wave case where the nodal lines give rise to a quadratic nonlinear contribution, we now find a nonlinear effect 
cubic in the applied field. Physically, this is quite transparent: the phase space volume available to the quasiparticle 
excitations increases as the cube of the field for point nodes, and as the square for line nodes. The nonlinear term 
anisotropically increases the magnetic field penetration because of quasiparticle occupation near the nodes, i.e., fewer 
Cooper pairs are participating in the current responsible for bulk flux exclusion. 

We can gain some insight into these results by examining the spatial de penden ce of the nonlinear part of the field as 
displayed in Fig. |l|. There the quantity H n i x , defined as the last term in ( 2.20a ) normalized to unity at its maximum, 



is plotted as a function of dimcnsionless distance D from the surface (D = Y sz — Y z ). The thickness of the sample 
is taken to be d » A z so that the behavior shown is that corresponding to a thick slab. The nonlinear field is 
constrained by the boundary conditions to vanish at the surface: the boundary condition implies an extremum for 
the nonlinear flow field at the surface and since H x oc dv± z /dY z , we see that H n i x must vanish there. It then increases 
rapidly reaching its maximum at about one half of a penetration depth and then decays exponentially inside the 
sample, as does the linear part. Thus arises the characteristic maximum of the nonlinear field seen in this Figure. 
The current is most easily obtained from H through Ampere's law for steady-state currents, which gives the result: 
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3z{Y z ) 



cH a sin tp 
Att\ z cosh(y 2S ) 



3mh ^ - 8 Uo 2 cosh(n„ 



fj(Y z ) 



(2.22a) 



c£f a COS Ip 



■ sinh^). 



Att\ x cosh^,,) 

The Y z dependence of the nonlinear current is contained in fj(Y z ) which is given by: 
fj(Y z ) = -3Y Z cosh(F z ) + (9/4) sinh(3F 2 ) - ( Cl + 6) sinh(n). 



(2.22b) 



(2.22c) 



In Eq. (2.22a), the first term on the right is the linear contribution, proportional to the applied field, while the second 
term is the nonlinear correction, modifying the total current. 

We now turn to the case where the OP is of the form ( |2.8| ). For si mplic ity, we consider tetragonal symmetry and 
quote the results for the orthorhombic case later. The current in Eq. ( [2.12 ) can be si mpli fied using that the magnetic 
field is applied in the x — z plane, so that v y is zero. The current is substituted in (2.9) to give, for the component 
containing the nonlinear term: 



d 2 v x 



0. 



(2.23) 



The method of solution is identical to that above with the only major difference being that the expansion parameter 
is now ol x — (vfx/^o)> linear rather than quadratic. The linear velocity field vq x is vq x = — c x sinh(y), while the 
nonlinear term is written as 



where 



vi x (Y x ) = (l/6)c x |c x | [cosh(2Y;) - 4cosh(Y x ) + 4c 2 sinh(Y x ) + 3] , 
e\ x H a cos ip 



ccosh(Y xs ) 



c 2 = tanh(y xs ) - sinh(Y xs ). 



(2.24) 



(2.25) 



The magnetic field is calculated again via the London equation. In this case the nonlinear part is (as in the d-wave 
case) proportional to H%, rather than to H% as was found for the point nodes. This follows again from phase-space 
arguments. After including the contribution from the purely linear component of v one finds: 



H Z (Y X ) = 



H a cos ip 
cosh(Y xs ) 



COsh^) + - 



1 / H a \cosip\ 



6 \H 0x cosh(Y Xi 



9h(Y x ) 



H X {Y Z 



H a simp 



cosh(F z ). 



cosh(y; s 

Here gn, which determines the nonlinear contribution to the magnetic field penetration, is given by: 

9h{Y x ) = -2sinh(2F x )+4sinh(y x ) -4c 2 cosh(y x ). 



(2.26a) 



(2.26b) 



(2.26c) 



It is again useful to plot the nonlinear component of the magnetic field. We consider H n i z , the last term in ( 2.26a ) 
normalized to its maximum value. Fig. ^ shows H n i z plotted versus dimensionless distance D from the surface 
(D = Y sx — Y x ). One sees again the rapid increase of the field near the surface of the sample, followed by the usual 
exponential decay. The plot is very similar to that in Fig. [I] for point nodes, but the field decays less rapidly into the 
sample. Within about three penetration depths, the nonlinear field is reduced to 20% of its maximum magnitude. 

The total current is composed of linear and nonlinear terms, as found from Ampere's law: 



3x(Y x ) 



Jz(Y z ) = - 



cH a cos ip 
4ir\ x cosh(Y xs ] 
cH a sin ip 



47rA z cosh(Y zs ) 



sinh(F x ) - 
sinh(Y z ). 



H a \ cos^l 
H 0x cosh(Y xs ) 



9J(Y X ) 



(2.27a) 
(2.27b) 



where 



9j(Y x ) = 4[cosh(2F x ) - cosh(Y x ) + c 2 sinh(F x )], 



(2.28) 



is the function determining the penetration of the nonlinear currents. 

We can now use our solutions for both OP's to derive expressions for the experimentally relevant quantities. 
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E. The transverse magnetic moment 



The expression for the magnetic moment in terms of the currents is:E3 

m=i /drrxj(v). (2.29) 
By making use of standard identities and the parity of v, Eq. ( 2. 29] ) can be rewrittenEl'El more conveniently as: 

(2.30) 



VH axz Acv zx 
m x . z = =F 



4~ 



2?re 



where A is the surface area of the plane along which the field is applied, V is the volume of the sample, and we have 
used that v is odd in z. The magnetic moment perpendicular to the applied field is given by m x cos ip — m z shiip. 
The linear terms of the velocity fields contribute to this quantity only if there is anisotropy in the penetration depth 
tensor. In that case, the linear term in the transverse magnetic moment, denoted by mj_, is: 



m ± = — AH a (X z - X x )sm2ip. 

47T 



(2.31) 



This term can be distinguished from the nonlinear contribution, m±, because of its different field and angular depen- 
dences. 



For an OP with point nodes, (2.7), m± is obtained from (2.30), (2.19) and (2.2C) as: 



(V.) = ^A\ z H a [(3 3 / 2 /32)/ sW ]/C 5 (F zs ). 



(2.32) 



This quantity is proportional to the cube of the applied field, rather than, as in the d-wave case,oE3€ZH23 to the 
square. This reflects the reduced phase space when the nodal regions are points rather than lines on the FS. Thus 
larger values of H a are very advantegeous provided that H a is kept below the field of first flux penetration, ff/i, so 
that the sample remains in the Meissner regime. The angular dependence of m±(ip) is contained in f s (ip) which is 
normalized to unity at its maximum and given by: 



f g (iP) = (16/3 3/2 )cos?/>sin 3 ?/>, 
while the dependence of m± on the material thickness d is given by the function ICs, 

K. s (Y zs ) = (l/2)sech 4 (r zs ) [3Y ZS - 2sinh(2F 2S ) + (1/4) sinh(4F 2S )] . 



(2.33) 



(2.34) 



The torque associated with m± is simply obtained by multiplying ( p. 32 ) by H a , therefore it has the same thickness 
and angular dependence. 

The function f s (ip) is displayed as the solid line in Fig. The transverse magnetic moment and torque in this case 
are maximal for the field direction corresponding to ip — tt/3, and vanish at directions corresponding to the nodes 
or antinodes of the OP. The it periodicity of m± matches that of the energy since the angular dependence of the 
quasiparticle energy arises solely from that of |A(0)| 2 . 

Since m± is an extensive quantity and it is often the case that larger samples can be made in film form rather than 
grown as free standing crystals, it is of considerable interest to examine the thickness dependence of the results, as 
given by ICs- The behavior of ICs is displayed in Fig. [|, where ICs is plotted (solid line) as a function of Y sz = dj (2A Z ). 
It is seen that ICs increases rapidly with Y sz , reaching 90% of its maximum value of unity when d ~ 5 A. If the sample 
is a thick slab, d 3> X z , then ICs — ► 1, so that m± is (for the same area A) maximal and independent of d. On the other 
hand, the decrease of ICs with thickness is substantial: for films where d — X z , ICs is /Cs(l/2) = .017, a reduction 
of over 80%. Such a decrease, however, may very well be compensated by a larger increase in A, compared to a free 
standing crystal. However, for extremely thin films, d -c X z , then ICs — l/4:0(xLJX z ) 5 , and despite the increase of Hfi 
in thin films, the amplitude of the signal would almost certainly be too smallnfl 



When the order parameter is of the form (2.8), the results obtained f rom the previous expressions (2.24) and (2.26) 
for a compound with tetragonal symmetry yield, when substituted in ( 2.30| ), the expression: 



m±(ip) 



^AX x H a Q|L) [|4/.r 2 )/Ac)}k: c iY,, 



(2.35) 



7 



We have introduced the ructions / C (V'), and Kc characterizing, respectively, the angular and the thickness dependences 
of the result. They are given by: 

f c (^P) = (3 3/2 /2)|cosV>|cosV>sinV', (2.36a) 
K c (Y xs ) = [sech(Y xs ) - l] 2 [1 + 2 sech(Y xs )} . (2.36b) 

The nonlinear transverse moment is now proportional to the square of the applied field, as in the d-wave case, because 
of the linear character of the nodal regions. The function f c (ip) is normalized to unity and it is plotted as the dashed 
line of Fig. |. It is seen there that the angular signature is different from that in the previous case: m± has now 
a maximum when ip — arctan( v / 2/2), although it is zero again for fields applied along the nodes and antinodcs 
(ip = 7r/2, 0) in the OP. As in the previous case, Kc is small for small thickness. If we are dealing with a thick slab, 
K c -> 1, but when d ~ A x , K c « .036. In the limit d < A x , Kc ^ 3/64(d/A :E ) 4 . As seen in Fig. ||, the overall 
characteristics of Kc (dashed line) are very similar to those of Ks, but Kc is larger in magnitude throughout. Both 

curves show a 50% drop in signal when d ~ 3X X . 

It is somewhat tedious but straightforward to generalize, starting from the form ( 2.12] ) , the calculation of the fields 



and currents for this order parameter to the case where there is penet ration depth anisotropy in the a — b plane. The 
result is that, when the sample is rotated about the b axis, Eq. ( |2.35 ) is simply modified to: 



miM = ^\ x H a (J^-) [(A/3 5 ^ 2 )f c miCc(Y xs )T x , (2.37) 

where T x = (l/5)(4 + S 2 ). Thus, the angular, field and thickness dependences remain the same, while only an overall 
anisotropy factor is needed. 



F. Penetration depth 

Measuring the field dependence of the penetration depth at low temperatures is another possible way of exploring 
the nonlinear Meissner effect. The reduction of the current via quasiparticle populationr-eesults in a lower superfluid 
density and hence a larger penetration depth. Indeed, this was the first quantity studiedta in this area, although it is 
only very recently that experimental measurementscil have been attempted for HTSC's. 

In the presence of nonlinear effects, several possible definitions of the penetration length which coincide in the linear 
limit give slightly different results. The appropriate definition depends on the experimental setup. We have calculated 
above the spatial current and field distributions, and these results can be used to obtain the nonlinear contributions 
to A for any definition. We briefly illustrate this here by computing the components of the penetration depth along 
the x and z directions via the definition: 



1 f dH x \ 1 1 / dH z 



A 2 (H a ) X z H a \ dY z J y=y A x (H a ) X x H a \dY x y ^ _j 



(2.38) 



where A^ is the zero field penetration depth along the i direction. We will assume that the field is applied along a 
symmetry direction for each OP and will not place any restri ctio ns on the thickness d. 



We consider first the order parameter with point nod es, (2.7 ). The penetration depth when the applied field is 
perpendicular to the z-axis, (ip = 7r/2) is obtained from ( [2.38 ) and ( 2.20| ): 



3 f H, s 2 



tanh(y zs ) - — — M C s }, (2.39) 



X z (H a ) A 2 I 32 \Hq 

where Cs = sech i (Y zs )[—4Y zs + sinh(4F zs )]. The most obvious difference between this and the results for d-wave is 
that the field correction is proportional to the square of the field, rather than to the field itself. This follows, once 
more, from phase space arg uments. For a thick slab, Cs — * 8, while in the very thin film limit, Cs ~ (4/3)(d/A z ) . 
For the OP given in (|2.8|) one similarly gets: 



Ax (H a 



1 = _L (tanh(y xs ) --( C c \ 
H a ) A x I. 3 \H 0x J J 



(2.40) 



where Cc = 1 — sech 3 (Y xs ), which goes to unity for a thick slab. For a very thin film superconductor, Cc — 2/3(d/\ x ) 2 . 
Here the penetration depth correction, 5X/X X is linear in the applied filed as in the d-wave case, as a result of the 
presence in both cases of nodal lines. The dependence on thickness in the very thin film limit is cubic in d for the 
point nodes and quadratic for line nodes. 
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III. DISCUSSION 



We have calculated the magnitude of the nonlinear electrodynamics effects as a function of field and angle. We 
believe that some remarks are now in order as to the feasibility of observing effects of the rough magnitude of those 
predicted. In making these remarks, we do not claim any experimental expertise in the relevant areas. We have in 
mind compounds with materials characteristics such as those of UBei3 or Sr2Ru04. In mentioning these compounds, 
we do not intend to propose that any of them belong to a specific pairing state. We merely wish to roughly estimate 
the level of the signal that would be predicted in the event that the material turned out to have a pairing state with 
a certain nodal structure. Our considerations can straightforwardly be extended to other materials. 

First- jOiir calculations have been performed, strictly speaking, in the low temperature limit. In practice, tfy 
meansE§L3 the temperature regime in the region T < T(H) < AqH/Hq, so that thermal excitations do not destroy 
the nonlinear Meissner effect. Assuming that, in order to maximize the signal, the applied field is close to Hfi, the 
field of first flux penetration, this still implies that the experiments must be performed at temperatures well below 
T c . This means, for the materials of interest, at dilution refrigerator temperatures. This is undeniably a disadvantage 
when compared wth the situation for HTSC's, but one that can be overcome by using torque magnetometryEJ or 
torsion oscillatorcJ techniques to measure the torque associated with the transverse moment. rXhese techniques can 
be adapted to use in conjunction with a dilution refrigerator and their sensitivity can surpassed that of the SQUID 
methods used for HTSC'al3a. These considerations pertain to the transverse moment and the associated torque. 
Several groups are planning, in the context of checking the very low temperature behavior of A in HTSC's for deviations 
from the linear power law behavior, measurements of A at dilutiQft,refrigerator temperatures. These techniques could 
be combined with the high resolution methods already employecH to measure the field dependence of A. 

We consider next the magnitude of the low T effect. The maximum amplitude of the transverse moment depends on 
the values of the penetration depth, the characteristic field Hoi, and the maximum field one can apply while remaining 
in the Meissner regime, which is Hn. Because m± is an extensive quantity, it also depends on the size (specifically 
the surface area) of the samples available. As an illustrative exercise, we have estimated a putatisiCp^iEnal for the 
transverse magnetic moment amplitude for a number of compounds by getting from the literatureBtj'Q 113 values of 
available crystal sizes and of the experimental parameters (such as penetration depths, and correlation lengths in the 
appropriate directions) that appear in our expressions. These values are subject to very considerable uncertainty and 
in most cases different references do not agree wi th ea ch othe r, bu t they are sufficient for the purposes of our exercise. 
By inserting them in the appropriate formulae, ( 2.32) ) , and ( 2.35| ), we obtain numerical values of the possible signal. 
The results are summarized in Table I, where we present our estimate for the maximum amplitude M± for several 
materials. M± is defined as the value of m-j for a thick slab at H a = H c \ and at the angle ^ for which m± is maximal. 
The critical field H c \ is calculated fromo (</>o/(47rA 2 )) ln(A/£) and is used as a conservative approximation 

to Hfi since H c \ is smaller (by a large factoiuia for YBCO) than -H/i. We have also included in the Table, for the 
purposes of comparison, one typical HTSC compound (YBCO), with the signal in that case computed from Ref. |28| 
for a d-wave state. For the other materials, we have assumed the OP in Eq. (2.7) for S^RuO/i, while for the listed 



heavy fermion materials and organic salt we have taken the OP of Eq. (2J5). It bears repeating that these choices 



are illustrative and do not imply any judgement on our part as to the likelihood of what the pairing state actually 
might be. Rather, our point is that the techniques in this paper can be implemented to infer the nodal structure of 
the pairing state. The results in Table I are expressed in physical units and also, for purposes of comparison, as ratios 
to the corresponding estimate for YBCO in a d-wave state. The numbers in the table are very encouraging: they 
are in all cases comparable to or larger thae, those for YBCO and always comfortably exceed the resolution of the 
experimental techniques discussed aboveESEJ For the penetration depth results, the situation is similarly favorable, 
since the changes in the penetration depth induced by a field close to Hfi are considerably larger than the lower limit 
(a few A resolution) already achieveded in YBCO. 

We have to consider also the limitations of this work and the presence of other phenomena, besides temperature 
excitations—tuat may reduce the signal. First, there is the question of impurities. As has been seen in the context 
of YBCOjOtJ good quality samples characterized by a transition temperature not appreciably degraded, and by 
the appropriate power law behavior of A with temperature, should exhibit a signal substantially of the magnitude 
calculated here for a clean system. The decrease in the nonlinear signal associated with nonlocal effects at lower fields 
can also complica.te the situationEJ. However, these effects are quite small for fields close to Hfi in the typical situation 
where this fielcO is considerably larger than the equilibrium H c \. In any case nonlocal effects are absent for several 
special crystal orientationsciH, which can then be chosen. There also important questions as to what the effect of 
using more realistic forms (still containing point or line nodes) of the order parameter would be, or of including in 
more detail the local value properties of the Fermi surface. All told, however, we believe that the estimates in the 
Table show that there is a sufficient cushion between the maximum value and the experimental resolution so that one 
can expect an observable signal. 
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In conclusion, we have calculated here the nonlinear signal arising from the presence of point or line nodes in some 
simple p-wave order parameters. We have shown that there is a likelihood that these effects will be observable in 
materials currently being studied. The results given can straightforwardly be extended, if and when the experimental 
situation warrants it, to the study of the low frequency responses, to more complicated or mixed order parameters, 
and a more general node spectroscopy procedure for p-wave materials can be performed as in Ref. E8[ 
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FIG. 1. Spati al beh avior of the nonlinear field when the OP has point nodes. The quantity ploted is H n i x , as given by 
the last term in (2.20a), normalized to unity at its maximum. This quantity is plotted as a function of D, the distance from 

D = Y sz — Y z . The sample has thickness d » A z . The nonlinear field increases rapidly until the 
and then decreases exponentially, as in the usual linear case. 



the surface in units of A 2 : D = Y sz — Y z 
distance from the surface is about half A 



FIG. 2. The nonlinear z comp onent of the magnetic field for a material with an OP of the form (2.5). The quantity plotted 
is H n i z , the last term in ( 2.26a| ), normalized to unity at its maximum. It is plotted versus dimensionless distance from the 
surface: D = Y sx — Y x . The material thickness is d » \ x . The behavior is qualitatively similar to that shown in the previous 
Figure, but the maximum occurs at a somewhat greater depth. 



FIG. 3. The angular dependence of the normalized transverse magnetic mom ent is plotted versus if) (the angle between H a 
and the z axis, in radians). The solid l ine is the function /s(V0 given in (2.33) for an OP of the form (2.7), while the dashed 
line is fc(ip) from (|2.3q) for the OP (2.8). Both functions are normalized to their maximum values. Their maxima are at 



tp = 7r/3 and if> = arctan(\/2/2) respectively and their periodicity is tv 



FIG. 4. The symbol K stands for the functions ICs (solid curve) and ICc (dashed curve), which characterize the thickness 
dependence of m±. They are plotted versus the dimensionless thickness Y s (Y s represents d/2X z for ICs and d/2X x for the 
ICc plot.) The exponential increase of the functions towards unity is rapid and the corresponding bulk regime arises when the 
material thickness d is about five penetration depths. 



TABLE I. Illustrative estimates of possible nonlinear signal for various materials of uncertain OP nodal structure. Values 
for d-wave YBCO are also given for comparison purposes. The magnetic field is in Gauss, and the maximum transverse 
magnetization amplitude, M±, as defined in the text, is given in (xl0~ 8 ) emu's. The quantities A and Ho are for the relevant 
directions (see text). 



Compound 


Area (mm 2 ) 


A (A) 


H 


H cl 
H 


M± 


M ± 
M | (YBCO) 


YBCO 


1.8 


1400 


7607 


.048 


2.1 


1 


UBei 3 


5 


11000 


204 


.032 


2.4 


1.1 


UPt 3 


150 


15000 


71 


.045 


66 


31.4 


Sr 2 Ru0 4 


25 


1940 


147 


.29 


22.5 


10.7 


(TMTSF) 2 C10 4 


1 


5000 


43 


.25 


2.7 


1.3 
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